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ABSTRACT 

The methods and techniques for the slitless spectroscopy software aXe, which was designed to 
reduce data from the various slitless spectroscopy modes of Hubble Space Telescope instruments, 
are described. aXe can treat slitless spectra from different instruments such as ACS, NICMOS 
and WFC3 through the use of a configuration file which contains all the instrument dependent 
parameters. The basis of the spectral extraction within aXe are the position, morphology and 
photometry of the objects on a companion direct image. Several aspects of slitless spectroscopy, 
such as the overlap of spectra, an extraction dependent on object shape and the provision of 
flat-field cubes, motivate a dedicated software package, and the solutions offered within aXe are 
discussed in detail. The effect of the mutual contamination of spectra can be quantitatively as- 
sessed in aXe, using spectral and morphological information from the companion direct image (s). 
A new method named 'aXedrizzle' for 2D rebinning and co-adding spectral data, taken with 
small shifts or dithers, is described. The extraction of slitless spectra with optimal weighting 
is outlined and the correction of spectra for detector fringing for the ACS CCD's is presented. 
Auxiliary software for simulating slitless data and for visualizing the results of an aXe extraction 
is outlined. 

Subject headings: Data Analysis and Techniques, Astronomical Techniques, Astrophysical Data 



1. Introduction 

A spectrum produced by a dispersing element, 
usually a prism or a grism, without a slit or 
slitlets has a number of characteristics that re- 
quire specialized data reduction tools. Slitless 
spectroscopy is primarily a survey tool, since the 
spectra of all objects within a given area defined 
by the telescope-instrumcnt-detector combination 
arc recorded. While this survey aspect brings ad- 
vantages, it also has some drawbacks. There can 
be no prior selection of target objects, and con- 
tamination, which is an overlapping of spectra, oc- 
curs frequently. The lack of defining slits implies 
that for extended objects the spectral resolution 
is controlled by the object size, more specifically 
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the component of the object size in the disper- 
sion direction. Moreover the shape of the objects, 
which in general are not aligned with the disper- 
sion direction, must be taken into account dur- 
ing the spectral extraction. Since in slitless spec- 
troscopy each pixel can receive radiation of any 
wavelength (within the sensitive range of the in- 
strument), special techniques for flat-fielding and 
fringing are required. 

Slitless spectroscopy has been employed for 
surveys from the ground such as for t he Ham- 



burg/ ESO objective-prism survey survey (jWisotzki et al. 
Il996t ). the UK Schmidt Telescope (Clowes et al. 
f980; Hazard et al . f986) and the Quasars near 
Quas ars survey (jWorseck. Wisotzki fc Selman 
120081 ). Early efforts to use this technique for 
the determination of stellar rad ial velocities are 
described in iFehrenbach ( 19691 ). However, the 
strong disadvantage of slitless spectroscopy from 
the ground is the high background, since the lack 
of a slit implies that each pixel receives the full 
transmission of the dispersing element. The back- 
ground can be reduced by choosing a high resolv- 
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Table 1 

The HST slitless spectroscopic modes to which aXe has been applied. 



Instrument/Camera Disperser Wav. Range [A] Dispersion [A/pixel] FOV 1 [arcsecond] 



ACS/WFC 


G800L 


5500- 


10500 




38.5 


202 x 202 


ACS/HRC 


G800L 


5500- 


10500 




23.5 


29 x 26 


ACS/HRC 


PR200L 


1700 - 


3900 


20 [ 


@2500A] 2 
31500A] 2 
@1500A] 2 


29 x 26 


ACS/SBC 


PR130L 


1250 - 


1800 


7[( 


35 x 31 


ACS/HRC 


PR110L 


1150 - 


1800 


10[ 


35 x 31 


NICMOS/NIC3 


G141 


11000- 


19000 




80.0 


51 x 51 


WFC3/UVIS 3 


G280 


2000 - 


4000 




13.0 


160 x 160 


WFC3/IR 3 


G102 


7800- 


10700 




25.0 


123 x 137 


WFC3/IR 3 


G141 


10500 - 


17000 




47.0 


123 x 137 



1 Field-of-View 

2 On account of the dispersion change with wavelength for prism spectra, the resolution is accompanied by 
the wavelength it refers to (see also Eqn. [3]) 

3 Prior to the WFC3 installation in Servicing Mission 4, aXe was used in the reduction of ground based 
calibration data. 



ing power disperser or adding a narrow band filter 
to the light path, such as for emission line studies 
on known lines (e.g. Salzer et al. 2000; Jangren et 
al. 2005). 

Deep slitless surveys from the ground, such as 
searches for high redshift galaxies at 3 < z < 7, 
are still valuable (see Kurk et al. 2004), how- 
ever space based observations are particularly ef- 
fective, since the background is many orders of 
magnitude lower than from earth and there are no 
strong and variable atmospheric absorption and 
emission components. For the Hubble Space Tele- 
scope (HST), the high spatial resolution (through 
lack of atmospheric seeing) brings a further gain, 
since the background contribution to the extracted 
spectrum is reduced. Although the background in 
slitless mode is larger by the ratio of the spectrum 
length to the slit width with respect to slit spec- 
troscopy, HST slitless spectroscopy can challenge 
the faint limits of gr ound-based telescopes (e.g. 
Malhotra et"al~ll2005h . 

Surveys with HST in slitless mode have in- 
cluded the NICMOS/Hubble Space Telescope 
Grism Parallel Survey of Ha emission line galaxies 
i n the redshift range 0.7 t o 1.9 with the G141 grism 
((McCarthy et al. I fl999h : the Space Telescope 



Imaging Spectrograph Parallel Survey (Gardner et 
al. 1998; Teplitz et al. 2003) of [O II] and [O III] 
emission line galaxies 0.43 < z < 1.7; and the 
Advanced Camera for Surveys GRAPES st udy of 
the Ultra-Deep Field (jPirzkal et al. 1 12004 with, 
for example, detectio n of Ly-a and Lyman- break 
galaxies to z ~ 6 ( Malhotra et al. 1 120051 ): the 
PEARS surv ey in the GOODS Nort h and GOODS 
South fields ( Malhotra et al. l2008f ) and searches 
for high-z supernovae to z ~ 1.6 (Riess et al. 2004; 
Strolger et al. 2004). 

After exploratory attempts to use standard 
spectral extraction tasks in public data reduc- 
tion packages (e.g., IRAF, IDL and MIDAS) for 
slitless spectroscopic data, it was decided to de- 
velop a dedicated extraction package. These stan- 
dard packages encapsulate an enormous amount of 
expertise for slit spectroscopy. However the lack 
of solutions for specific problems in slitless spec- 
troscopy together with the demand for the pro- 
cessing of hundreds to thousands of spectra per 
image mandated a dedicated software package. 

This software package, called aXe, was origi- 
nally developed for Advanced Camera for Surveys 
(ACS) grism and prism spectroscopy (see Table [T] 
and Pasquali et al. 2006). However by design aXe 



2 



was intended to be instrument-independent, and 
this aim was achieved by placing all the instru- 
ment specific parameters into a single data file 
called a configuration file. In addition to ACS, 
aXe has subsequently been successfully appli ed to 
NICMOS slitless data (|Freudling et al. 112008? ) and 
data obtained during the ground calibration cam- 
paigns of Wide Field Camera 3 (WFC3, Bond et 
al. 2007, Kuntschner et al. 2008b; see also Table 
UK- 

This paper is devoted to a description of the 
aXe spectral extraction package. aXe does not 
work like a data reduction pipeline, which has a 
static set of products as output. Depending on 
the data and the scientific questions to be an- 
swered, the user can reduce the data following dif- 
ferent recipes and methods, and the subsequent 
Sections illustrate all the methods and techniques 
implemented in the current version aXe-1.6. Sec- 
tion [5] is devoted to the extraction of spectra from 
an individual slitless image. Section [3] explains 
the various techniques to treat the contamination 
problem, followed by the introduction of a 2D- 
resampling technique which aims at co-adding in- 
dividual spectra of the same object in Section 0] 
The aXe implementation of optimal weighting is 
illustrated in Section [5l followed by a discussion 
of the effects of detector fringing (Section [6]). The 
structure of aXe and related software packages for 
simulation and visualization are presented in Sec- 
tion [7J The paper closes with a summary and 
options for further refinement of the package. 

2. Extracting spectra from single slitless 
images 

This section details the aXe method to extract 
spectra from individual slitless images. aXe re- 
quires associated direct images for the slitless im- 
ages to enable data processing. 

2.1. Slitless image direct image relation 

In slitless spectroscopy the location of spectra 
on the detector is defined by the position of the ob- 
ject itself within the Field-of-View (FOV). There- 
fore, in order to recover the absolute wavelength 
scale, the position of the object must be known in 
advance in the reference frame of the slitless im- 
age. Although not a fundamental limitation, the 
presence of a direct image from which the position 



of the object can be determined is usually taken 
as a pre-requisite for the reduction of slitless spec- 
tra. In instrumental configurations with a grism 
as disperser, the th order spectrum, when avail- 
able, could be used to provide the zero point of the 
wavelength scale. However, due to the prism on 
which the grating is ruled, the th order appears in 
practice slightly dispersed which does not favour 
an accurate position determination. 

In order to establish for each object a reference 
position on the slitless image, aXe requires, for 
each slitless image (Figure [TJ panel a), a corre- 
sponding direct image of the same sky area (Fig- 
ure [TJ panel b). The pointing information as given 
in the world coordinate system information in the 
headers of the slitless and direct images allows to 
transform object positions from the direct images 
to the slitless images. The reference positions are 
then used as the origins of local coordinate systems 
to describe, for each object, its spectral properties. 

2.2. The aXe configuration and calibration 
files 

aXe has been applied to the various HST slit- 
less spectroscopic modes listed in Table [1] which 
all have different properties, e.g. different descrip- 
tions of the trace (the location of the center of 
gravity in the spatial direction) and dispersion. In 
order to be able to reduce these distinct instru- 
mental modes with aXe, all information for a spe- 
cific mode is stored in a single configuration file. 
An aXe configuration file contains: 

• the data description (pointers to the science, 
error and data quality FITS image exten- 
sions); 

• the location and extent of the different spec- 
tral order(s); 

• the trace description for the spectral or- 
der (s); 

• the dispersion solution for the spectral or- 
der (s); 

• the names of the calibration files (flat-field 
and sensitivity files) to be used. 

The trace description in aXe is given as a poly- 
nomial: 



Ay(Ax) = ao + ai * Ax + 02 * Ax 2 



(1) 
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Fig. 1.— Cutouts (1200 x 700 pixel) of an ACS/WFC G800L slitless image and the associated F850LP 
direct image are shown on panels a and b, respectively to illustrate the input data for aXe. Bad pixels and 
cosmic ray hits have been removed. Panels c and d show intermediate results from a spectral extraction with 
aXe. The local background image (see Sect. 12.6. 2|) is shown in panel c and a model of the slitless image, 
which is used to estimate the contamination (see Sec. 13. 2|) in panel d. Since the extraction which produced 
images c and d was limited to bright objects (juab < 22.0), some fainter sources are not masked out in the 
background image and are not part of the model image in panel c. The +l st (A), th (B), +2 nd (C), +3 rd 
(D), —1 st (E) and — 2 nd (F) order spectra for two objects 1 and 2 are marked in the model image (panel d). 



with (Ax, Ay) — (x — x re f,y — y r ef), the rela- 
tive distance of the image coordinates (x, y) from 
the reference position of an object on the slitless 
image (x re f,y re f). Grism dispersion solutions are 
described as polynomials: 



A(0 = /o + *i*^2*/ 2 



(2) 



with the distance along the trace Since prisms 
display quite rapid variation of the spectral dis- 
persion with wavelength, the functional form for 
the dispersion solution is better described as: 



A(Z) = h + 



k 



Is 



(l-h) (l-l a ) 2 (l-l ) 3 



(3) 



Around the singularity at / = Iq the resolution de- 
creases rapidly, causing a pile-up in instrumental 
configurations such as ACS/HRC/PR200L which 



still have significant s ensitivity at the start wave- 



2006) 



lengt h of the pile-up (jLarsen. Walsh fc Kiimmel 



All quantities do, cti, a%, . . . , Iq, l\, I2, ■ ■ ■ in Eqns. 
[T][3]can be given as a 2D field dependent polyno- 
mial in order to take variations of the trace de- 
scription and dispersion solution as a function of 
object position into account. This 2D field depen- 
dence is caused by field effects and the different 
geometrical distortion between the direct and the 
dispersed images. The quantity a given as a 2D 
polynomial of order 2 is, for example: 



a + 

ai * x + ct2 * y + 

a3 * x 2 + 04 * x * y + as * y 2 



(4) 




the extraction geometry (the direction of 
constant wavelength and slit width, see Sect. 
I2.4p and the extraction area are determined; 



Fig. 2. — The projection of pixel pj onto the 
set of spectral bins [. . . , 6»_i, bi, ■ ■ ■]■ For the 
given trace and extraction direction (gray, thin 
lines), pj contributes with the weight (= area) 
ajj-uctj^ajj+i to the bins bi, b i+ i. 



6. 



where x and y are the reference object position on 
the slitless image. 

A thorough description of the aXe configuration 

files and the various quantities declared therein 

can b e found in the aXe Manual (jKummel. Kuntschner fc Walsh 
2008). As calibration files aXe needs sensitivity 
files to transform the extracted spectra from de- 
tector units (electrons / s) to flux units (fx) and 
special flat-fields, which are introduced and de- 
scribed in Sect. 12.51 The HST instrument science 
teams conduct dedicated observing programmes to 
derive calibration files, trace description and dis- 
persion solution in the various slitless modes (see 
e.g. Pasquali et al. 2006; Kuntschner, Kiimmel 
& Walsh 2008a). aXe configuration and cal- 
ibration files for all supported instruments are 
available for download on the ST-ECF webpages 
(http: / /www. stecf.org/instruments/) . 

2.3. Spectral extraction 

The spectral extraction in aXe aims at avoiding 
multiple rebinning of the pixels and processes the 
data from detector to the ID spectra in a direct 
manner. This is achieved by tying the spectral 
information to the pixels. In detail, the aXe ex- 
traction of a single object spectrum from a slitless 
image involves the following steps: 



3. for each pixel pj in the extraction area, the 
distance to the trace in the cross-dispersion 
direction, the distance along the trace and 
the effective wavelength associated to the 
pixel are computed; 

4. each pixel pj is flat-fielded according to its 
associated wavelength (see Sect. I2.5|) : 

5. a set of spectral bins [6q, b\, b%, • • • , b n ] is cre- 
ated. Each bin bi covers the wavelength 
range 

AA = [(A(i) - A(i - l))/2, (A(i + 1) - A(i))/2] (5) 

where X(i) is the grism or prism dispersion 
solutions given in Eqns. [2] and [3] The se- 
quence of integer numbers i is chosen such 
that the entire sensitive wavelength range of 
he instrument is covered; 



the area of each pixel pj is projected onto 
the set of spectral bins bi in order to de- 
termine the fractional contributions a^i. As 
can be seen in Figure O the quantities aj,% 
are usually zero for all except a few bins. 
The extraction or cross-dispersion direction 
is optimized for each object (see Sect. I2.4[) 
and is, as shown in Fig. [2J not necessarily 
perpendicular to the trace; 



7. the spectrum value of the spectral bins with 



non-zero 



contributions 



are re-calculatcd 



using weighted summation 



bi 
4 



* S Pj + Wbi * Sb t 



(6) 
(7) 



where w' b ., s' b . and w^, are the new and 
old values for the spectrum value and weight 
for bin i, respectively and s Pj is the spectral 
flux in pixel pj . 



1. on the basis of the reference position, the 
local trace and dispersion solutions are de- 
termined; 



After the extraction the spectral bins for each 
object, orde red in wavelength , are written to a 
FITS table (jPonz et al. Ill994h . These resulting 
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Fig. 3. — An illustration of the parameters for 
defining the extraction geometry. The object mor- 
phology is described as a 2D Gaussian with major 
and minor half axes a and b and the angle O, here 
defined with respect to the dispersion or trace di- 
rection. The quantities p a and pb are the projec- 
tion of the major and minor half axis onto a plane 
perpendicular to the trace. 

spectra contain for each spectral bin the wave- 
length A, the extracted spectrum in [electrons / s], 
the error in [electrons / s], the flux in f\, the flux 
error and the total weight. 

The wavelength values, although derived from 
spectral bins bi which have associated wavelength 
ranges as described above, are not just assigned 
the mean value of the spectral bins X(i) (see Eqn. 
[5|), but are computed using weighted summation 
as in Eqns. [5] and [7| (with X Pj , the effective wave- 
length associated to pixel j, instead of s Pj ). 

The error for each spectral bin is computed by 
propagating the initial error for each detector pixel 
through all aXe reduction steps using basic error 
propagation. If not explicitly given in an image 
extension of the slitless image, aXe derives an error 
image on the fly, using a typical detector noise 
model (readout noise and photon shot noise). 

The total weight of a spectral bin is the co- 
added area of its pixel contributions. Since aXe 
does not use pixels which are marked as defective, 
the weights for a spectral bin can vary. For HST 
images, such bad pixels (e.g. dead pixels or pixels 
affected by cosmic rays) are indicated by the value 
of the data quality image extension. 



2.4. Extraction geometry 

Every object which produces a slitless spec- 
trum acts as its own virtual slit (see Appendix 
in Freudling et al. 2008 for a detailed discussion), 
which in general is not aligned with the detector 
grid or the spatial direction of the disperser. The 
shape of this virtual slit must be taken into ac- 
count if a ID spectrum is to be extracted from the 
slitless image or a rectified 2D spectrum (analo- 
gous to a long slit spectrum, see Sect. [4j is to be 
generated. Since there is no slit or mask which de- 
fines or limits the extraction geometry, the 'slit ori- 
entation' and the 'slit length' translate in slitless 
spectroscopy to the extraction direction (the direc- 
tion of equal wavelength) and the extraction width 
(the number of pixels on both sides of the trace 
which are co-added), respectively. These quanti- 
ties must be set for each object during the spectral 
extraction. 

aXe offers several methods to define the extrac- 
tion geometry and thus to extract ID spectra from 
a slitless image. Most of the methods link the ex- 
traction geometry to the object shape in order to 
optimize the extraction. Figure [3] illustrates the 
quantities involved in defining the extraction ge- 
ometry for an arbitrary object. In a very basic 
model an object is described as a 2D Gaussian 
with major axis size a, minor axis size b and the 
orientation of the major axis 8, which is defined 
with respect to the trace here. The projections of 
a and b onto the plane perpendicular to the dis- 
persion or trace direction are named p a and pb, re- 
spectively. For this object shape, aXe offers spec- 
tral extraction with: 

1. a fixed extraction width w = const, given by 
the user and an extraction direction perpen- 
dicular to the dispersion direction; 

2. a variable extraction width scaled as w = 
k * MAX(p a ,pb) (k given by the user) and 
an extraction direction perpendicular to the 
dispersion direction; 

3. a variable extraction width scaled as w = 
k*a (k given by the user) and an extraction 
direction parallel to a; 

4. as in [3J but switching to w — k * b and an 
extraction direction parallel to b if pb > p a - 
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Fig. 4. — The four coefficient images of the 
HRC/G800L flat-field cube. In the region (top 
left) obscured by the coronographic finger the co- 
efficients result in a flat-field value 1.0 for all wave- 
lengths. 

Depending on the science case the user may pre- 
fer a specific extraction method, e.g ([1} if the pri- 
mary targets are point-like or ((J) in a typical sur- 
vey situation with galaxies of differing sizes and 
random orientation. 

2.5. Flat-fielding slitless spectroscopic data 

Flat-fielding is an essential data reduction step 
in creating scientific quality images of uniform sen- 
sitivity. In direct imaging, flat-fields for each filter 
are constructed using homogeneously illuminated 
images, and a single flat-field correction yields 
equal sensitivity. These direct imaging flat-fields 
are valid for radiation in the restricted passband 
defined by the combination of optical elements. A 
slitless spectrum, however, can occur anywhere on 
the detector, and the flat-field correction has to be 
performed for each pixel as a function of incident 
wavelength; this leads to the use of a flat-field cube 
whereby the wavelength behaviour of the flat-field 
is described for each pixel by the 3 rd dimension in 
addition to the two spatial dimensions of the de- 
tector. To characterize the wavelength dependent 



behaviour, a series of direct image flat-fields taken 
at various central wavelengths is used. 

For ACS, we decided to combine the wavelength 
information in all direct image flat-fields by fitting 
a polynomial to the series of wavelength depen- 
dent flat-field values for each pixels. The flat-field 
cube then contains the polynomial coefficients in 
the third dimension. This method is applicable 
for the ACS slitless modes since the variation of 
the flat-field as a functio n of wavelength is smooth 
(jBohlin fc Hartig||2002| ). 

There are several reasons why the flat-field cube 
derived from direct image flats only does not pro- 
vide a perfect correction. The large scale flat-field 
characteristic of the grisms and prisms is likely 
somewhat different from that of the direct im- 
age flats used to generate the flat-field data cube. 
Moreover the direct image flat-fields produce a flat 
image for a homogeneously illuminated image of 
the sky, but the plate scales and hence the effec- 
tive pixel size vary signi ficantly from one corne r of 
the detector to another (jCox &: Lindler 1120021 ) . 

In order to correct for these additional, field 
dependent effects it is necessary to apply an em- 
pirical large scale correction to the flat-field cube, 
which is determined from observations of a flux 
standard star at various detector positions. The 
final flat-field cube allows the application of a sin- 
gle sensitivity curve for the conversion to absolute 
flux units which is independent of the object po- 
sitions. 

A detailed description on flat-fielding in slit- 
less spectroscopy and a discussion on the use 
of the available direct image flats is given in 
Walsh fc Pirzkal I (|2005l ). Figure 2] shows a repre- 
sentation of the ACS HRC/G800L flat-field cube. 
Its sub-panels show the four layers with the poly- 
nomial coefficients as indicated in the labels. The 
flat-field value for the pixel (i, j), which is exposed 
to light of wavelength A is computed as 

FF(i,j,x) = a (i,j) + a 1 (i,j)x+ (8) 
a2(i,j)x 2 + a 3 (i,j)x 3 

where cik(i,j) is the coefficient from the k th layer 
of the flat-field cube, x the normalized wavelength, 



x — (A — \min)/(^min — Amax) 



(9) 



and Xmin and Xmax are the minimum and max- 
imum wavelengths of the direct image flat-fields 
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Fig. 5. — The flat-field values as a function of 
wavelength for two ACS/WFC pixels (CCD #1) 
derived from direct image flat-fields. The coeffi- 
cients of the polynomial fits (shown as lines) to 
the data are stored in the flat-field cube. 



used in constructing the slitless flat-field cube, re- 
spectively. 

Figure [5] illustrates the flat-field as a function of 
A for two selected pixels on the ACS/WFC CCD 
#1. The symbols show the values from the direct 
image flats taken through the ACS filters for a 
2x2 binned region, with the rms on the mean as 
errors. The lines show the polynomial fits to the 
the data and thus the wavelength dependent flat- 
field correction which is applied to the data values 
measured in these two pixels. 

2.6. Background determination 

aXe has two different methods to handle the 
background emission from the sky. The first 
method incorporates a global background sub- 
traction before starting the extraction process, 
the second uses local sky background estimates. 
Further enhancements to the basic methods de- 
scribed here, e.g. applying a 2D correction to 
the master sky when th e background shape differs 
(jMalhotra et al. 1120081 ), are possible. 



2.6.1. Global background 

The method of removing a global sky back- 
ground uses a 2D master sky image which is au- 
tomatically scaled to the background level of each 
individual slitless image. When determining the 



scaling factor, pixels on and close to the spectral 
traces are masked out in order to prevent an over- 
estimate of the background level. The spectral 
extraction is then conducted after subtracting the 
scaled master background from the slitless images. 
Master sky images can be generated by e.g. com- 
bining multiple slitless images with a rejection al- 
gorithm which excludes pixels dominated by ob- 
ject spectra. 

2.6.2. Local background 

The local background estimate uses the back- 
ground pixel values close to the spectral traces 
to determine individual background levels. For 
each pixel which is part of an object spectrum, 
a background value is determined by interpolat- 
ing between the user-defined number of adjacent 
background pixels. Figure Q] shows the cutout of 
a local background image in panel c. Interpolated 
background pixels are visible at the position of 
the spectra, the background pixels used for the in- 
terpolation remain unchanged, and pixels neither 
part of a spectrum nor used as background appear 
white. 

The local background subtraction can even be 
applied in addition to the global background. Such 
a procedure is useful in slitless modes such as NIC- 
MOS G141, which display temporal variations of 
the background (see Freudling et al. 2008). 

2.7. Discussion of the extraction method 

The aXe approach to extract spectra is differ- 
ent from the long slit or multi-slit approaches. The 
method of equipping the pixels with spectral infor- 
mation (e.g. wavelength, trace distance, see Sect. 
12.3ft and then 'filling' them into initially empty 
spectral bins has several advantages. Firstly, the 
re-use of an individual pixel in several spectra, 
which is necessary in regions of spectral overlap, 
is straight forward. Several wavelength values 
Ao, Ai, A2, . . . are given to the pixel pi when co- 
added to the spectra of the sources 0, 1,2,.. .. Sec- 
ondly, the virtual slits defined by the objects them- 
selves can be flexibly addressed via the various ge- 
ometrical methods of attributing wavelength val- 
ues to the pixels (Sect. 12. 4| . Hence aXe is well 
equipped to adjust to the variety of object mor- 
phologies encountered in typical survey scenarios. 
Thirdly, the techniques used in the aXe spectral 
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extraction are geometrical descriptions (e.g. for 
the trace and wavelength solutions or attribut- 
ing wavelength values to the pixels) and weighted 
summation (when co-adding pixels to the spectra), 
which both are rather simple and mathematically 
well described. 

3. Contamination 

In slitless spectroscopy an overlap of spectra 
occurs very frequently. This overlap is generally 
referred to as 'contamination', as it complicates 
the task of extracting a clean spectrum of a given 
single object. There exist several types of con- 
tamination. Firstly, objects close in the spatial 
direction can cause contamination. Such contami- 
nation can also occur in spectroscopy with slits or 
masks. Secondly, objects that are well-separated 
on the sky can still contaminate each other in 
the dispersion direction. In spectroscopy with 
masks, such contamination is avoided by a care- 
ful object selection or, equivalently, mask design. 
Thirdly, for grism dispersers (e.g. ACS/G800L, 
NICMOS/G141) many spectral orders (0 th , 1 st , 
2 nd , —I s * . . . ) are typically visible in the FOV of 
the detector, and contamination can occur across 
different spectral orders at distances of several 
hundred pixels. 

Contamination is impossible to avoid in typ- 
ical slitless spectroscopy observations at opti- 
cal/NIR wavelengths, such as the HST slitless 
modes ACS/WFC/G800L, ACS/HRC/G800L, 
NICMOS/G141. Consequently, the data reduc- 
tion software must be able to collect and process 
contamination information in order to deliver data 
that can be interpreted scientifically. 

aX e handles contamination in tw o different 
ways (jKiimmel. Larsen fc Walsh I l2005h . Both use 
predictions based on the knowledge of the source 
positions. The so called geometrical contamina- 
tion scheme simply records which spectral bins 
are affected by contamination. The quantitative 
contamination scheme estimates for each spectral 
bin the contribution from neighbouring sources. 

3.1. Geometrical contamination 

For geometrical contamination the areas that 
are covered by the various spectral orders of all 
objects are marked. Then for every pixel in each 
spectrum the information on how many other 



spectra fall on this pixel is processed through the 
extraction. The final result is the number of other 
spectra that fall on the pixels which were com- 
bined to a spectral bin. Geometrical contamina- 
tion is computationally cheap and allows the selec- 
tion of regions of spectra which are not occupied 
and therefore contaminated by other spectra. 

3.2. Quantitative contamination 

Quantitative contamination aims at computing 
directly for every spectral bin the contaminating 
flux from other sources. This is done by mod- 
elling the dispersed contribution of every object 
to the slitless image. Based on this model infor- 
mation the contaminating flux for each pixel is 
recorded and processed through the spectral ex- 
traction. The result of quantitative contamina- 
tion is, for each extracted spectrum, an associated 
spectrum of the contaminating flux. 

aXe generates a slitless image model re-using 
the configuration and sensitivity files from the ex- 
traction, which guarantees a symmetry between 
the extraction and the modelling in quantitative 
contamination. The morphological and spectral 
object properties are derived from the direct im- 
ages, and aXe has two different methods to employ 
these properties in an emission model, which are 
discussed in the next Sections. 

An even better refinement and better estimates 
would be expected from an iterative approach 
which uses the direct image information initially, 
but then the extracted spectra for more refined 
contamination estimates. This procedure, which 
is computationally intensive and does not neces- 
sarily converge, is currently not implemented. 

3.2.1. Gaussian emission model 

The Gaussian emission model approximates the 
object morphology with 2D Gaussians, whose ma- 
jor and minor axes and position angle are already 
available in the spectral extraction (see Sect. 12.41 
and Fig. [3]). As spectral energy distribution this 
emission model uses one or, if available, several 
AB-magnitudcs at different wavelengths measured 
on the associated direct images. Beyond the wave- 
length range covered by the AB-magnitudes the 
object 'spectra' or spectral energy distributions 
are extended as flat in f\. Figure [TJi shows the 
slitless model generated for quantitative contami- 
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3.3. Contamination example and discus- 
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Fig. 6.— Extracted ACS /WFC/G800L spectra 
(solid lines) and quantitative contamination esti- 
mates (dotted lines, Gaussian emission model) for 
two spectra extracted with aXe. The contamina- 
tion estimate in the lower spectrum indicates that 
both emission lines results from an overlap with 
the spectrum of other sources. 



nation using the Gaussian emission model. 

3.2.2. Fluxcube emission 

The Fluxcube emission model derives both the 
object morphologies and the spectral information 
directly from the direct images. For this purpose a 
so called fluxcube file is generated for each slitless 
image. Such a fluxcube file is a multi-extension 
FITS image which contains at least one corre- 
sponding direct image cutout, converted to fx, and 
a segmentation image. The segmentation image, 
essential to avoid self-contamination, marks the 
pixels contributing to each object. The values of 
the pixels on which an object in the associated 
source catalogue was detected contain the object 
index. A seg mentation image is easily produced 
by SExtractor (jBertin fc Arnouts" 19961 ) (parame- 
ter ' CHECK I MAGE_TYPE SEGMENTATION'). As for the 
Gaussian emission model, the fluxcube emission 
model extends the object 'spectra' as flat in fx- 



Figure [5] illustrates how quantitative contami- 
nation aids in the interpretation of extracted spec- 
tra. The two panels show source spectra (solid 
lines) and quantitative contamination information 
(dotted lines, Gaussian emission model) for two 
objects extracted from HST ACS/WFC/G800L 
slitless data. The upper source (named ID4245 
internally) has a strong emission line at ~ 7300A, 
and the lower (ID1785) two emission features at 
~ 6600 and ~ 8600A. The geometrical contam- 
ination scheme raises the contamination flag for 
all spectral bins in both spectra. The quanti- 
tative contamination scheme is able to estimate 
for the upper spectrum a negligible contamina- 
tion from other sources and for the lower spec- 
trum at least a strong contributions from neigh- 
bouring objects to both emission features. The 
two spectra in Fig. [B] clearly demonstrate the need 
for quantitative contamination information. This 
is especially essential for deep slit less surveys such 
as GRAPES (jPirzkal et al. Il2004h . where contami- 
nation by th orders often resembles emission lines 
(as in ID1785) and a reliable distinction is essen- 
tial. 

The basic data necessary to extract slitless 
spectra with aXe (slitless images plus associated 
direct images) already allows the usage of quanti- 
tative contamination with little extra effort. A 
limitation of quantitative contamination is the 
emission model, and particularly its spectral part. 
The crude approximation of the object spectra 
with often only a single data point in flux space 
can still lead to large differences between the con- 
tamination estimate and the actual contamination 
in the spectra. E.g. an examination of the spectra 
TD1 785 in Fig. [6] reveals that the contamination 
estimate is ~ 2 to 3 times lower than the excess 
in flux at ~ 6600A. A better wavelength coverage 
with more direct images yields better contamina- 
tion estimates. 

4. 2D-resampling of spectra 

As in direct imaging, a typical slitless spectro- 
scopic dataset taken with HST usually consists of 
several images of the same regi on on the sky with 



only small shifts or dithers (jKoekemoer et al. 
2002) between them. The basic aXe extraction 
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discussed in Sect. [2] is capable of reducing all ob- 
jects from every slitless image. However the data 
for each object must be co-added to obtain deep 
spectra with high signal-to-noise ratio. 

4.1. Co-adding individual spectra 

Co-added object spectra can be derived by com- 
bining the ID spectra extracted from the individ- 
ual slitless images. This approach has some dis- 
advantages. The data is rebinned twice, once for 
the spectral extraction and once for combining the 
spectra. The procedure requires a complex weight- 
ing scheme to take into account masked pixels (e.g. 
cosmic ray hits) and different exposure times for 
the spectral bins. Information in the spatial direc- 
tion, e.g. the exact location of an emission feature 
within an extended source, is lost when combining 
ID spectra. Problem detection and error tracking 
is more difficult in ID than in 2D. 

4.2. Co-adding spectra with aXedrizzle 

With aXedrizzle (Kummcl ct al. 2004; Kiimmcl 
et al. 2005) we have developed and implemented 
in aXe a method of combining the individual, 2D 
spectra into a deep 2D slitless spectrum, which 
allows subsequent ID extraction from the deep 
image. For 2D-resampling and co-addition of the 
i ndividual 2D spec t ra we use the drizzle software 



(|Fruchter fc Hook I l2002h . which is a standard 
method of combining HST direct images. With 
aXedrizzle the re-gridding to a uniform wavelength 
scale and a cross-dispersion direction orthogonal 
to the dispersion direction is achieved in a single 
step. The weighting of different exposure times 
per pixel and cosmic-ray affected pixels are cor- 
rectly handled through the drizzle weights. The 
combined 2D spectra can be visually inspected for 
problem detection. They reveal fainter features 
in the spatial direction not detectable in single 
images (see e.g Rhoads et al. 2004) and can be 
used for finding emission line regions in larger, 
extended objects. 

4.3. The aXedrizzle technique 

As an example, Figure[7jillustrates the aXedriz- 
zle procedure for ACS/HRC G800L data from the 
Hubble Ultra Deep Field (HUDF; proposal ID 
9803, P.I: Beckwith; 

http://www.stsci.edu/hst/udf) Parallels observa- 
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Fig. 7. — An overview of aXedrizzle steps and 
products. The flat-fielded object pixels are re- 
assembled to a cutout image (a). Drizzling this 
image leads to an individual 2D spectrum with 
zero trace angle, linear dispersion and perpendic- 
ular extraction direction (b). The individual 2D 
spectra are co-added to a deep, 2D slitless image 
(c). 

tions for one object. The upper panel in Fig. [7J 
shows the cutout image for one object in one grism 
image, assembled from the flat-fielded pixels. In 
the 112 slitless images of the entire dataset, the ob- 
ject is located at different positions, which leads to 
small variations in trace angle, dispersion solution 
and image distortions. For each of the 112 obser- 
vations of this object individual drizzle coefficients 
are determined which, via the drizzle process, lead 
to individual 2D grism images, such as Fig. [7b, 
with zero trace angle, a linear, pre-defined disper- 
sion solution and an extraction direction perpen- 
dicular to the trace. aXedrizzle allows the user to 
select a specific interpolation kernel for the driz- 
zle routine in the configuration file. The final ID 
spectrum is extracted from the co-added 2D slit- 
less stamp image Fig. [7h, which has a combined 
exposure time of \2Aksec. 

To illustrate the advantages of aXedrizzle, Fig- 
ure[8]shows the aXedrizzle combined, deep 2D slit- 
less images corresponding to the spectra in Fig. [6] 
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Fig. 8. — The aXedrizzled, deep 2D slitless images 
for the sources shown in Fig. [6] For ID1785 (lower 
panel) the contaminating th order spectra from 
other sources are clearly visible. 

The re-sampled image confirms the intrinsic origin 
of the emission line feature for the source ID4245 
in the upper panel and even indicates the presence 
of additional emission lines. In the lower panel, 
it is obvious that the object ID1785 itself has a 
smooth spectrum. Both 'emission features' arise 
solely from two bright, offset spots, which are iden- 
tified as contaminating zeroth order spectra from 
two different objects. 

Besides the main science information, aXedriz- 
zle also processes other data such as the errors, 
contamination (see Figs. [5] and [5} and weights (for 
optimal weighting, see Sect. [5]). All data is stored 
in the layers of the deep 2D slitless stamp, which 
is a multi-extension FITS image, and used in the 
ID extraction. As in the original drizzle method, 
the weight image is composed as one image exten- 
sion which reflects the total exposure time for each 
pixel. 

The 2D resampling and co-adding in aXedriz- 
zle introduces correlated errors between ne i ghbor - 



ing pixels. However in Freudling et al. ( 2008T) 



it was shown that the independently processed 
noise from aXedrizzle is in good agreement with 
empirical noise estim ates derived from DEfLSNR 



(jStoehr et al. Il2008f ). This confirms not only the 



accuracy of the errors from aXedrizzle but also jus- 
tifies their further use, e.g. in optimal weighting 
(Sect.0. 

Local background estimates are processed by 
aXedrizzle in the same way as the object data, 
and for each deep 2D slitless stamp a 'deep' 2D 
background stamp image is generated from the 
local background estimates on the individual im- 
ages. As in the processing of individual slitless 



images (Sect. I2.6p . the background subtraction is 
done during the ID extraction of the spectra. 

Only the 2D combining of the first order spec- 
tra is currently implemented in aXedrizzle. In 
all grism modes supported by aXe, the sensitiv- 
ity contrast between the first and all higher orders 
is very large, and the focus on the first order does 
not imply an important loss of scientific content 
which the slitless data contains. 

4.4. Data with different roll angles 

Apart from dithering, HST slitless data can also 
be observed with different roll angles, which trans- 
lates to a different orientation of the traces on the 
sky. Certainly the simple method of co-adding ID 
spectra (Sect. |4TT|) . but also aXedrizzle, is capable 
of combining such data. As well as the spectra, 
the individual contamination patterns, which usu- 
ally change for significant roll angle differences, 
are combined. In general, this leads to combined 
spectra with larger and more complex contami- 
nation. Also care must be taken when combin- 
ing data of extended objects, since the object size 
in the spectral direction, and thus the spectral 
resolution, varies for significant roll angle differ- 
ences. In addition, more complex objects, such 
as late-type galaxies with HII regions, can display 
very different spectra at different roll angles. For 
point sources, however, it is possible to combine 
the contamination-free parts of the data only and 
to derive contamination free spectra. 

5. Optimal weighting in aXe 

The opti mal weighting tech nique, originally in - 
troduced bv lHorneT(|l986l ) and Robertson ( 1986I ). 
aims to enhance the signal-to-noise ratio of spec- 
tra by attributing lower weights to pixels which 
contain only a small fraction of the object flux. 
The basic equation of the spectral extraction us- 
ing optimal weights is (see e.g. Rodriguez-Pascual 
et al. 1999): 

er(x,A) 2 

where x is the coordinate in spatial direction, A 
the coordinate in the spectral direction, f(x, A) the 
data value at pixel (x, A), b(x, A) the background 
value at pixel (x, A), cr(x, A) the noise value at pixel 
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(x, A), p(x, A) the extraction profile at pixel (x, A) 
and /(A) the extracted data value at A. 

In the original descriptions of optimal weight- 
ing, the extraction profile p(x, A) is computed from 
the object spectrum itself. For several reasons 
this approach is not generally applicable in slit- 
less spectroscopy. The signal-to-noise ratio of 
the sources is often too low to determine an in- 
dividual extraction profile. The contamination 
phenomenon does not permit an automatic and 
reliable generation of extraction profiles for all 
sources. An iterative approach on the sometimes 
hundreds or even thousands of spectra on an ACS 
slitless image would require too much computing 
time. Also averaging the profile in the spectral 
direction does not give representative results for 
special object classes, e.g. sources with emission 
lines. 

In aXe we have developed an optimal weighting 
which is robust and practical for all sources. To 
compute extraction profiles, the optimal weighting 
in aXe uses the models for the dispersed objects, 
which were introduced in Sect. 13.21 as the basis of 
the slitless models and hence the quantitative con- 
tamination technique. The source-specific models 
computed from the direct image (using either the 
Gaussian or the Fluxcube emission model) are a 
good basis for calculating the quantity p(x, A) in 
Eqn.riOl 

The 2D models are also used as an input to 
calculate the pixel errors a(x, A) according to the 
typical detector noise model: 



a (x. A) = y mod(x, A) + b(x, A) + rdnoise 2 (11) 

where mod(x, A) and rdnoise are the 2D model 
value at pixel (x, A) and the detector readout 
noise, respectively. 

6. Treatment of Detector Fringing 

Optical interference in a CCD occurs between 
incident light and that back-reflected from the thin 
multi-layer structure. In a back-illuminated CCD, 
if the detection layer becomes transparent, the 
contribution of the reflected light can become a 
large fraction of the incident power. The result is 
a modulation of the signal by fringes whose ampli- 
tude depends on the optical properties of the lay- 
ers. The level of fringing and its behaviour with 



wavelength is dependent on the number of layers of 
the CCD, their thicknesses and refractive indices. 

In imaging and slit spectroscopy fringe maps 
are used to correct for the detector fringing. In 
the case of slitless spectroscopy, parallel to the case 
for flat-fielding, each pixel can receive light of any 
wavelength over the passband of the instrument, 
and the fringe map must be replaced by a cube. In 
practice maps of the individual layer thicknesses, 
together with their optical properties, are used to- 
gether with a simple geometrical optics model to 
predict the fringe amplitude for all p ixels for any 
wavelength. Malumuth et al. I ( 2003h applied this 



technique to the STIS SITe CCD and an analo- 
gous procedure was developed for the A CS WFC 



and HRC detectors ( Walsh et al. 1 12003) 



From a set of monochromatic flat-fields over 
the range of significant fringing, the Fresnel equa- 
tions can be solved in terms of the thickness of the 
CCD layers, adopting tabulated refractive indices 
for the layer materials (jWalsh et al. 1120031 ). In or- 



der to correct for the fringing at a given pixel in 
a slitless spectrum, the wavelength extent of each 
pixel (given by the position of the dispersing ob- 
ject and the dispersion solution), the CCD thick- 
nesses for that pixel and the respective refractive 
indices are used to compute the observed fringing 
of that pixel. 

The correction of the fringe amplitude affecting 
a slitless spectru m has been incorporated into the 



aXe software and iKuntschner. Kiimmel fc W alsh 



(|2008al ) present the application to ACS slitless 
spectra. Besides the CCD structure, an impor- 
tant input for the fringe calculation is the form 
of the throughput function for the CCD pixels. 
The pixel throughput function is defined as the 
relative distribution of light as a function of wave- 
length falling onto a given CCD pixel. In general 
this is the convolution of the intrinsic source spec- 
trum with the line spread function (LSF) of the 
instrument. The determination of the throughput 
function and the computation of the fringing am- 
plitudes would thus require the a priori knowledge 
of the source spectrum, which in general is not 
available. However for continuum sources and the 
theoretical case of pure emission line objects (an 
object with a spectrum that consists only of a sin- 
gle, narrow emission line) the throughput function 
can be estimated. 

For continuum sources the LSF dominates the 
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form of the throughput function everywhere and 
the influence of the slowly varying source spectra 
can be neglected. Using a Gaussian with 100A full 
width at half maximum height (FWHM) as LSF 
(similar to the spectral resolution, see Kuntschncr, 
Kummel k Walsh 2008a) the ACS WFC requires 
a maximum fringe correction of 0.1%. This is neg- 
ligible in comparison with other sources of error 
(photon noise, flux calibration errors). 

For pure emission line sources the throughput 
function is solely determined by the intrinsic ob- 
ject spectrum, thus amplitudes of 12% as mea- 



surcd in the monochromator data ( Walsh et al 



erate on datasets, which usually consist of several 
slitless images and their associated direct images 
taken in the same area of the sky. 

Some preparatory steps are required before 
spectra can be extracted. These steps mostly deal 
with generating input source lists from the direct 
images in the dataset. The direct images are usu- 
ally co-added to enhance the signal-to-noise ratio 
and to remove pixel defects and cosmics. For 
HST data, this step is gener ally performed us- 



ing the MultiDrizzle software (IKoekemoer et al 



120061 ) . An object detection algor ithm, e.g. SEx 



tractor ( Bertin k Arnouts 19961 ). is run on the 



2003) and even higher are possible. While such 
variations could be detected, in practice the in- 
trinsic line broadening in the spectrum dilutes 
the fringe signal. As an example, line flux varia- 
tions of < 5% were detected among a set of Wolf- 
Rayet star observations with ACS WFC /G8000L 
(|Kuntschner. Kiimmel k Walsh I [2008ah . If the 
emission line source is extended, then mixing of 
wavelengths in the resulting slitless spectrum will 
even further dilute the detected fringe amplitude. 

7. The aXe software package 

The aXe extraction software is available for 
download. In this section we describe the oper- 
ational system of version 1.6, together with two 
related software packages which are used to simu- 
late and visualize slitless data. 

7.1. aXe extraction 

aXe is distributed to be installed and loaded as 
a local IRAF/PyRAF package 

(http: / /www. stecf.org/ software/slitless software/axe/) . 
Moreover it is embedded into the STSDAS (http: / /www 
software package which is distributed by the Space 
Telescope Science Institute (STScI). The aXe soft- 
ware package is written in ANSI C and Py thon It 



1999). 


GSL 


(Mink 
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forms. 



(GoughJ 120061) and WCSTOOLS 



The aXe software is split up into several tasks, 
which perform all the necessary steps to extract 
spectra from slitless images. There exist two dif- 
ferent classes of aXe tasks. Low Level Tasks oper- 
ate on individual images, and all input and output 
refers to a particular image. High Level Tasks op- 



co-added direct image. Besides detecting the ob- 
jects, SExtractor is also capable to deliver all nec- 
essary object parameters such as the positions, 
Gaussian shapes (for defining the extraction ge- 
ometry) and source brightness (for quantitative 
contamination). Using SExtractor is not manda- 
tory, however the source lists for aXe must be 
formatted in the same way as a SExtractor object 
list. The positions of the sources determined on 
the co-added direct image are first projected back 
into the coordinate system of each associated di- 
rect image, and aXe provides a dedicated task to 
do this step in a typical scenario of a direct im- 
age co-added with MultiDrizzle and a SExtractor 
source list. Then the object positions are pro- 
jected into the coordinated system of the slitless 
images as part of the aXe extraction. 

7.2. Simulations with aXe 

The quantitative contamination method intro- 
duced in Section 13.21 is based on modelling the 
slitless spectra on input images, either as 2D 
Gaussians (see Fig. [TJi) or object images. Thus, 
st^p^Sjikkng this modeling capability to a full sim- 
ulation tool was a natural development. Rather 
than adding simulations to aXe, an additional 
PyRAF package called aXeSIM 

(http://www.stecf.org/software/slitlesssoftware/axesim/; 
see Kiimmel, Kuntschner k Walsh 2007), which 
uses a number of the capabilities of aXe already 
described, was built as a tool for simulating slitless 
spectroscopic images plus their associated direct 
images. 

aXeSIM was designed both to assist users plan- 
ning slitless observations and to aid in the analysis 
of slitless data, such as for determining detection 
limits by adding and extracting model spectra on 
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Fig. 9. — aXeSIM simulations of bright objects 
(m FW ow < 23.0) in the NICMOS HUDF field 
for the HST/WFC3 G141 grism and F160W filter 
in the upper and lower panel, respectively. The 
WFC3/IR FOV is 123" x 137". 



real images. For HST, aXeSIM is a valuable tool 
in the Phase I and II proposal preparation pro- 
cess, and its usage has already started in Cycle 
17. A 2D impression of the layout of a target field 
is provided, which allows assessment of potential 
problems arising from crowding or spectral over- 
lap. The sensitivity and resolution offered by the 
various slitless modes can be explored in 2D, such 
as for a variety of object morphologies. The op- 
timal pointing and roll angle can be chosen for a 
specific target field. The characterization of the 
instrument used by aXeSIM through a configura- 
tion file is identical to that used by the extraction 
package aXe, thus closing the loop between the 



simulation and subsequent extraction. 

In the most basic form an object is simulated in 
the Gaussian emission model (see Section l3.2.ip by 
a Gaussian shape and one AB-magnitude, trans- 
formed and extended as a spectrum flat in f\. For 
simulating more realistic objects, the user can pro- 
vide 2D image templates for real object shapes, 
build more complex spectral energy distributions 
by specifying magnitude values at different wave- 
lengths or provide high resolution spectra as tem- 
plates, which are shifted in redshift and scaled in 
flux to user-provided values. The basic input for 
every object is collected in a SExtractor-like text 
table. 

As a typical application of aXeSIM, Figure 
[9] shows a noise-free simulation of the bright 
(m FW0W < 23.0) objects in the NICMOS Hub- 
ble Ultra Deep Field (Thompson et al. 2005) 
for the HST/WFC3 G141 grism and F160W fil- 
ter in the upper and lower panels, respectively. 
The spectral energy distribution of the objects 
is defined via their AB-magnitudes in the fil- 
ters ACS / WFC /F850LP, NICMOS3/F110W and 
NICMOS3/F160W, which all were provided in the 
original source catalogue. 

7.3. aXe2web visualization 

A deep slitless image (e.g. from ACS/WFC) 
can contain detectable spectra of hundreds to 
over a thousand objects. To visually inspect 
them in a fast manner we developed aXe2web 
(http: / /www. stecf.org/ software/slitless-software/axe / axe2web.php) , 
a tool which produces browsable web pages for fast 
and disc erning examination of m any hundreds of 
spectra (jWalsh & Kiimmel |[2004 . 

This additional task to the aXe package takes 
the aXe output files and generates linked web 
pages in html format for all sources. Each object 
produces a line in a table which lists the reference 
number, magnitude in the direct image filter, the 
X and Y position of the direct object, its Right 
Ascension and Declination, a cut-out image show- 
ing the direct object, the spectrum stamp image 
showing the 2D spectrum, a ID extracted spec- 
trum in counts and the same in flux units, over- 
plotted with the 'spectrum' of the contaminating 
sources. Figure [10] shows an example of two ob- 
jects extracted from an ACS/WFC/G800L grism 
dataset as presented by aXe2web. 
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Fig. 10. — An example the visualization of aXe 
products with aXe2web. The columns of the html- 
table show some basic object information (se- 
quence number, ID number, brightness, position) 
followed by a direct image cutout and a 2D slit- 
less image. The extracted spectra are plotted in 
[electron/ s] and in physical units (/a), overplotted 
with the contamination estimate. 



8. Conclusion and outlook 

aXe is not the only software for extracting HST 
slitless spectra. For ACS/WF C, an independent 



metho d has been d escribed in iDrozdovskv et al. 
(|2005l) . MULTISPEC dApellaniz Il2005l) and the ex- 



tractio n method described in lChen. Lanzetta fc Pascarelle 
( 1999h concentrate on STIS slitless spectroscopy. 
Quickly after the first release in 2002, aXe devel- 
oped into the standard reduction package for all 
ACS slitless modes. The methods implemented 
in aXe specifically for slitless spectroscopy (e.g. 
flat-fielding, 2D-resampling with aXedrizzle) have 
proven to produce science-ready and well cali- 
brated data for a number of large pr ogrammes 
such as GRAPES dPirzkal et al. I l200l. PEARS 
([Malhotra et al. 1120081 ) and searches for high-z su- 
pernovae (Riess et al. 2004; Strolger et al. 2004). 

By generating instrument-specific configuration 
and calibration files, aXe can also be applied to 
slitless spectroscopic data from other HST instru- 
ments (e.g. NICMOS) and even from ground 
based instruments such as the Wide Field Im- 
ager (WFI) mounted on the ESO 2.2-m-telescope. 
With small modifications, aXe has even been used 



to successfully reduce data taken with FORS2 
MXU masks (see Kuntschner et al. 2005 and 
Kummel et al. 2006). aXeSIM is being used to 
simulate data for future space missions such as 
SPACE (Cimatti et al. 2008; now renamed to EU- 
CLID) and SNAP (http://snap.lbl.gov/) that plan 
to use slitless spectroscopy. 

The ability within aXe to automatically extract 
spectra from slitless data was important in the 
NICMOS Hubble Lega cy Archive (HLA) project 
([Freudling et al. 1 120081 ) and aXe served as the 
central component of the completely automatic 
Pipeline for Hubble Legacy Archive Grism data 
(PHLAG) (Kummel et al. 2007; Kummel et al. 
2008). Due to the modular setup of aXe, ad- 
ditional tasks which handle some features spe- 
cific to NICMOS could comfortably be added in 
order to achieve optimal results leading to an 
archive of several thousand near-infrared spectra 
(http: / /www. stecf.org/ archive/hla/). 

Besides instrument-specific adjustments also 
the core of the aXe extraction software is contin- 
ually being improved and further developed, with 
releases of new versions at an approximately an- 
nual rate. New features and developments that 
are currently being implemented for aXe (versions 
1.7 and 1.8) and aXeSIM (version 1.2) include 
spectral extraction based on a slit geometry opti- 
mized for slitless spectroscopy and flux conversion 
which compensates for the size of the direct ob- 
ject (developed in the NICMOS HLA project, see 
Freudling et al. 2008) ; the identification and exclu- 
sion of cosmics and hot pixels within aXedrizzle , 
similar to MultiDrizzle ( Koekemoer et al. 20061 ) 



in direct imaging; the integration of aXeSIM into 
STSDAS. 

The application of aXe to HST data is going to 
continue with the installation of the Wide Field 
Camera 3 (WFC3, Bond et al. 2007) with its three 
slitless spectroscopic modes (see Tab. during 
Servicing Mission 4. As for ACS, aXe is the rec- 
ommended tool for reducing WFC3 slitless spec- 
troscopic images. From the ground calibration 
data collected during the Thermal Vacuum test- 
ing, initial aXe configuration and calibration files 
are already available (Kuntschner et al. 2008b and 
Kuntschner et al. 2008c) that were, for example, 
used in generating the aXeSIM simulations shown 
in Fig. [9l Improved calibration will result from on- 
orbit data, enabling aXe to extract spectra from 
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all slitless modes of WFC3. 

Based on observations made with the NASA/ESA 
Hubble Space Telescope, obtained from the data 
archive at the Space Telescope - European Coor- 
dinating Facility. 

We would like to thank S.S. Larsen, W. Freudling, 
and M. Dcmlcitncr for their respective contribu- 
tions to aXe. We are grateful to our anonymous 
referee for the numerous suggestions, which helped 
to significantly improve this paper. 
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